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Abstract 

The temperature dependence of the Casimir effect for the radiation field confined between 
two conducting plates is analysed. The Casimir energy is shown to decline exponentially with 
temperature while the Casimir entropy which is defined in the text is shown to approach a limit 
which depends only on the geometry of the constraining plates. The results are discussed in terms 
of the relation between the Bose distribution function and the equipartition theorem - a relation 
based on a study by Einstein and Stern circa 1913. 



1 Introduction 



The Casimir effect involves the difference between the energy of a field subject to some constraining 
boundaries and the energy of the field with such constraints removed. Although Casimir Q introduced 
this consideration long ago, in 1948, the study of this effect enjoys currently considerable popularity 
[||, |3|, [|, [||. Indeed this journal included very recenty an exposition on the effect ||. The experimental 
status of the effect was tenuous until 1997 with the report of vindication of the theoretical prediction 
to an accuracy of a few percent The reason for this continuous interest is not hard to find: in 
this effect we have a direct demonstration of the existence of the so called "vacuum fluctuations" or 
"zero point energy" . Thus the original effect was a prediction that two parallel conducting plates will 
experience mutual attraction in vacuum (at zero temperature, T = 0), caused by modification of the 
allowed modes due to the presence of the plates. Remarkably the attractive force is independent of 
the coupling of the electromagnetic (EM) field to matter ( viz., the electronic charge, e) and is propor- 
tional to the velocity of light, c and Planck's constant, h. Today, of course, we understand that every 
quantized field exhibits fluctuations even at its lowest (vacuum) state and these fluctuations sample 
the allowed modes. Furthermore, whereas the original effect related to modified boundary conditions, 
more recent works (e.g., Schwinger as cited in ||) led to studies of bulk (volume) Casimir effect (where 
we consider the forces on a sample because it is different from its environment, e.g., through having a 
different dielectric constant) and to the dynamical Casimir effect where time is involved. 

Another direction for Casimir effect studies is its extensions to finite temperatures. Again the sim- 
plest case deals with modified modes of the EM field due to the presence of conductors as boundaries 
but now each mode is thermally occupied in addition to it being sampled through vacuum fluctuations. 
This also leads to an effective force between the plates. Finite temperature studies were carried out 
by several authors H. Yet, to our knowledge, the classical limit of the effect, even for the EM field, 
has not been clarified todate. This paper addresses this problem. 

The paper is organized as follows. The next section, section 2, contains a review - hopefully self 
contained - of the Casimir effect problem in its simplest form, viz., radiation field confined between 
two conducting plates. This section includes the exact solution for the Casimir energy, free energy 
and entropy |l(J. These are displayed in a figure which gives the behavior of these quantities for all 
temperatures. The succeeding section, Section 3, includes the central point of the paper, viz., an expo- 
sition of the classical limit and a study of its implications. Section 4 is devoted to the demonstration 
that the classical solution is "robust" - by this we mean that a naive high temperature expansion for 
the Casimir energy doesn't have any nonvanishing corrections to the classical value. In the conclu- 
sions which are given in Section 5 we review the old (cf. |2j) argument for the existence of zero point 
fluctuations as viewed from our vantage point. 
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2 The Casimir Effect 



Evaluation and definition of the Casimir energy at zero temperature due to vacuum fluctuations for 
the case under study here is given in several texts and reviews (e.g., H g 0, Hi). Our presentation 
therefore, although aspiring for self containment, is somewhat sketchy. We consider the radiation field 
confined between two conducting plates. The size of the plates edge is L. The first plate is placed 
at z = in the XY plane, and the second at z = d parallel to the XY plane. L » d. (In fact 
we are interested in L — > oo while d remains finite.) The energy, i.e. the expectation value of the 
Hamiltonian, tied down in zero point fluctuations in the mode k is, in obvious notation (k includes 
the polarization and the zero in the argument relates to T = 0), 
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The total energy density of the EM field with the conducting plates as boundaries is (fen is the 
magnitude of the wave vector parallel to the plates 0] ) 
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The energy density, in dimensionless units, e(d,0), is given by 
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The equation above defines D whose dimension is energy density. The energy density due to vacuum 
fluctuations of the radiation field in an arbitrarily large volume, V (e.g., V = L 3 ) - which serves as 
the reference, unconstrained, system is 
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e(oo,0) = J xdx J dm\J x 2 + m 2 , 
The Casimir energy density, at T = 0, in dimensionless units, is given by 

e c (0) = e(d,0) -e(oo,0). 



(6) 



Both e(d,0) and e(oo, 0) diverge. e c is commonly [Q, ||, |Tl| evaluated by a physically justifiable 
regularization technique. Thus a wave vector dependent function, r(k/k c ), is introduced into the 
above integrals such that for k >> k c r — > while r — * 1 for k « k c thereby rendering the integrals 
convergent. A simple choice is (a = l/k c ), 

r(ak) = exp [—ak] . 

a is allowed to go to zero at the end of the calculations - for all the terms together - thus the sum is 
"regularized" . The details of the calculations will not be given here (cf . |TJ] ) . This issue is further 
discussed in section 4. The result for our case is 

sc(0) = -^C(4), (7) 
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where 
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The force per unit area, T '/L 2 , between the plates is now calculable from the regularized Casimir 
energy, Eq. 7 , yielding, after reverting to physical dimensionality, 
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thereby yielding the well known Casimir force n]. 



E(d,0) -£7(oo,0)] 



240d 4 ' 



(8) 



The finite temperature problem is quite similar. Now the zero point energy is supplemented by 
the thermal energy. Thus the (average) energy tied down in the mode labeled by k is 



E k {T) = (l/2)fiLo k + 
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Correspondingly, 
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Returning to our dimensionless units we may write the energy density as a sum of a zero temperature 
part plus a temperature dependent part: 



e(d,T) = e(d,0)+u'(d,T), 
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The corresponding expressions for the unconstrained system are 

e(oo, T) = e(oo, 0) + «'(oo, T), 

with 

oo 

u'(oo,T) = / dm f{m). 



(13) 
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Evaluating the sum in Eq. 12 via the Poisson summation formula [11, 10 1 gives (fi = 2irm) 

OC -j^ rp 

e c (T) = e(d, T) - e(oo, T) = -4t 3 ^ -coth^csch 2 ^), t = vr— . 

m=l ^ c 

For t << 1, e c — > e c (0). 



(15) 



The finite temperature Casimir energy, e c (i), i.e. the expectation value of the Hamiltonian of the 
constrained system with the that of the unconstrained subtracted from it, is displayed in Fig.l. The 
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relevant free energy was calculated by several authors ||. We outline the procedure as follows: The 
partition function (for one mode) is 



Z k = > exp [-f3huj k {n + 1/2)] = — - — -, 

^ 1 - exp [-PhoJkl 



(16) 



and hence the free energy for the mode labeled by k is 



Ft = -k B TlnZ k = -huk + k B T]n(l - exp [-(3huk\), (17) 



Hence the expression for the free energy density of the constrained system is, 
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A corresponding equation holds for the unconstrained system thereby leading to the dimensionless 
expression for the Casimir free energy density, (ft c , given by 

^ = DMt), (19) 
where D is defined by Eq. 3, and (detailed derivation is given in the appendix), 

cj)Jt) = _2t V — \co\h(tn) + (^)csch 2 (£//)] . (20) 

m=l ^ 

The Casimir free energy, <p c {t), is displayed in Fig. 1. For t « 1, (f> c (t) — > e c (0), as it should. 

At this juncture it is natural to consider Casimir's Entropy [jicfl . This entropy, cr c (t), in our 
dimensionless units is defined by 

4>c{t) = e c {t) - ta c {t). (21) 

— a c (t) is displayed in Fig. 1 (note the negative sign). At high temperatures, t » 1, s c {t) falls off 
exponentially with t, cr c {t) approaches a constant value (independent of t and, of course, of K) while 
4> c (t) becomes proportional to t. Reexpressing these results in standard physical dimensions we get, 
at this limit, t » 1, 

E c ^0, F C ^-TS C , Sc ^^(^). (22) 

Note that the entropy is proportional to the area of the plates scaled by d 2 (d is the plates' separation). 
We thus conclude that in the classical limit (high temperatures) the Casimir force is purely entropic. 
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3 The Classical Limit - Discussion 

An intuitive understanding of the vanishing of the Casimir energy for T » T c ( cf. Fig. 1) may 
be gained via the following reasoning. Loosely speaking || the number of normal modes (per unit 
volume) in our confined system is unchanged upon changing d (in the case under study). (A proof of 
this, based on regularization, is given in the next section.) Thus moving the walls adiabatically leads 
to shifts in the levels - not to appearance (or disappearance) of levels (modes). This (coupled with 
the classically valid equipartitioning of energy, which means that each harmonic oscillator- like mode 
holds ksT amount of energy) implies that the energy density, in this classical limit, is unchanged upon 
changing d. i.e., in the classical limit, defined by the validity of the Rayleigh Jeans (RJ) law or the 



equipartioning of the energy [14|, the Casimir energy is nil. This, since it is defined as the difference in 
the energy (density) between the constrained (d « L) and the unconstrained {d ~ L) cases and we 
have just argued that there exists (intuitive) one to one correspondence between the levels regardless 
of the size of d. Thus we have, in conformity with the result of the exact calculations as exhibited in 
Fig. 1, that, 

lim e c (d,T) = 0. (23) 



T 

Tr 



In discussions of black body radiation, Kirchhoff 's law (that the ratio of emissivity to absorptivity of 
all bodies is a universal function of the wavelengh, A and the temperature, T) is used § to infer that 
U, the total electromagnetic (EM) energy in a cavity at thermal equilibrium is a function of T only, 
U = U(T) with U the thermodynamic internal energy. We see from our discussion above that this is 
true only for T » T c where T c is a characteristic temperature of the cavity (cf. Eq. 13). (e.g., at 
T = it is, strictly speaking, never true.) For this reason we refer to Eq. 23 as "Kirchhoff 's theorem". 



We remark that the RJ result can be readily obtained within classical physics [14] and hence our 
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definition of classical limit as the one where the equipartition theorem holds is a reasonable one. 



The entropy, S, as a function of the energy, E, of a cavity (i.e., a constrained system), is given by 
(e.g., El) 

S(E) = -k B ln £(£), (24) 

where X(-E) is the number of states with energies less than or equal to E. In evaluating the Casimir 
entropy in the classical limit (which is the difference in the entropies, at high energies, of the entropies 
of the constrained and the "free" systems) we expect the dominant contribution, the "volume term" 
[||, to cancel. This is because at these energies the dominant contribution is from the short wavelength 
modes (A ~ 1/T). These are insensitive to the boundaries when the dimension of the cavity much ex- 



ceeds these wavelengths [15|. The contribution to the Casimir entropy, then, is essentially determined 
by the long wavelengths which are geometry dependent, i.e., which relate to the "shape" of the cavity 
and is independent of temperature. It follows then that the Casimir free energy is proportional to the 
temperature. This reasoning is born out by the explicit calculational results. Hence we see that this, 
viz., the free energy being proportional to T, is not at all related to the RJ law - the latter relates to 
the energy. 



We are now in a position to interpret the zero point energy (zpe) as a contribution necessary to 
assure that at high temperatures , the energy, U, is a function of T only as, indeed, was noted long 
ago (1913) by Einstein and Stern (||, p. 2): without it the energy will depend on the boundaries. Al- 
ternatively, if we assume the validity of "Kirchhoff 's theorem" at high temperatures, we may deduce 
the zero temperature Casimir energy as follows. 



Let the energy per mode, k, of the radiation field be written as (T c refers to the constraints, if 
present) 

u(k, T; T c ) = u{k, 0; T c ) + u'(k, T; T c ) 



where v! is the energy held in the mode without the zpe. For an allowed mode, k: 
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The Casimir energy is then, 



e c (T;T c ) = Y,u(k,T;T c ) -J2<k,T). 



(25) 



(26) 



(27) 



Use of Eq. 25 ( i.e., separating the thermal energy from the vacuum's) and Kirchhoff 's theorem ( Eq. 
23) implies, 



lim 

T/T c ^oo 
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. k k 
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with e c (0) being the Casimir energy at zero temperature. Hence 



lim 

T/T c ^oo 



u'(k, T; T c ) - ^2 u'{k, T) 
. k k 



-e c (0). 



(28) 



(29) 



Now both terms in the square bracket are readily calculable - their temperature dependence assures 
convergence - and thus yield the Casimir energy at T = ( e c (0) ) without recourse to a (perhaps) 
objectionable regularization scheme. 
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4 Robustness of the Classical Limit 



In this section we evaluate the Casimir energy through a power series expansion in t (t = vr^-). This 
scheme requires regularization for each term in the expansion. The result is that, within such an 
expansion, Kirchhoff 's theorem is exact, i.e. the Casimir energy vanishes to all orders in this high 
temperature expansion. This is interpreted as implying that the classical equipartition theorem is ro- 
bust - the classical approximation, once taken, is exact to within power series corrections. We hasten 
to add that this is yet another example of incorrect handling of infinities which are further discussed 
at the end of this section. 



Let us return to the expression for the Casimir energy density at finite temperature, Eqs. 10 and 



15, which, in our dimensionless units are summarized by (x m = V x 2 + m 2 ) 
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We now adjunct to each of the above integrals the "standard" [11, 12 1 cutoff function, 

/ = exp [-ax m ] , 

assuring thereby the convergence of the integrals. ( We are interested in the a — » limit.) We note that 
such regularization scheme is justifiable on physical grounds as follows [11, 12, The conductivity 
is, in general, a function of frequency. Indeed, all metallic conductors are effectively transparent 
to radiation with wave lengths comparable to the interatomic spacing, a. Hence the integrals (and 
summations) over wavenumbers, k, are limited to k < k c ~ 1/a, i.e., for wavenumbers \k\ > k c the 
plates do not provide a boundary and hence the Casimir energy for k > k c evidently vanishes. 
Assuming the validity of the "standard" regularization scheme (each integrand is multiplied by a cutoff 
function assuring the convergence of the integrals) we may expand flq| , p. 35, 
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£?2n are Bernouli numbers. 



The evaluation of £ c (t) reduces to the evaluation of terms like [1C] 
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The cutoff parameter, a, is set equal to zero at the end of the calculations; the result is 
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Noting that fij, p. 1076, 



y 
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we get g p = for p even, and for p odd, 
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for p = n — 3. 



Returning to our expression, Eq. 30, we see that upon substituting Eq. 31 only terms with p even 
occur, i.e., s c (T) = 0. This confirms the exponential decay of the Casimir energy with temperature 
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that is implied by Eq. 15, and depicted in Fig. 1. (The case of p = will be recognized as the sum 
over the constrained number of modes, per unit volume, with the unconstrained number deducted 
therefrom. That the result is nil constitues a proof, based on a particular regularization scheme, of 
our assersion (Section 3) that the number of modes, per unit volume, is unchanged upon varying the 
plates' separation, d.) 

This "robustness" of the classical (erroneous) solution is worthy of note and it implies that the 
correct expression can't be obtained from the wrong one by analytic (power series) means. 



5 Concluding Remarks 

This paper gives the classical - here the high temperature - limit of the Casimir effect. The case 
studied was the simplest and historically the first considered - viz. the radiation field confined between 
two large (in the limit - infinite) parallel conducting plates separated by a relatively short distance, 
d. We argued that in the classical limit (defined to be temperatures such that the Rayleigh Jeans, 
i.e. energy equipartition, theorem holds) "Kirchhoff's theorem" is valid - i.e. the energy density 
of the radiation field is a function of the temperature only and this implies that, in this limit, the 
Casimir energy vanishes. We showed that the zero point energy is required to assure the validity of 
"Kirchhoff's theorem" in the classical limit. Alternatively, assuming the validity of the theorem allows 
the evaluation of the Casimir energy at T = without recourse to any regularization scheme. We 
noted, following reference H, that these results were anticipated by Einstein and Stern in 1913, prior 
to the formulation of quantum mechanics and quantum field theory. Thus these authors noted that 
the high temperature expansion of what we now term Bose distribution function is (fihu « 1) 

t^s 1 - k B T - hco + 0(1/T). 
exp(pTkj) — 1 2 

Thus we have here a temperature independent term which contributes to the total energy. Its can- 
cellation, i.e., the validity of "Kirchhoff's theorem" in our presentation, requires a positive zero point 
energy, -\-\lvjj. Thus the removal of the zero point energy by considering a "normally" ordered Hamil- 
tonian does not eliminate the need for zero point energy. It (the zero point energy) is seen, ala the 
Einstein and Stern result as discussed above, to be required for the correct classical limit or, in our 
terminology, for the validity of "Kirchhoff's theorem". This stems here from the form of the Bose 
distribution rather than the Hamiltonian. 



The evaluation of the (Casimir) free energy in the classical limit led to showing that the Casimir 
entropy, which is defined in a natural way in the text, is (in this limit) temperature independent 
constant and reflects finite volume corrections in statistical physics. Thereby the proportionality of 
the free energy to the temperature, T that is well known @ is seen not to be related to the Rayleigh 
Jeans limit [17] but rather is a consequence of the perhaps interesting result that the Casimir force in 
the classical limit is purely entropic - in fact geometric. We have demonstrated that the classical limit 
results, i.e., at temperatures high enough to warrant the applicability of the Rayleigh Jean's law are 
"robust" in the sense that (naively) it leads to vanishing corrections. 
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Appendix: Derivation of Equation (20). 
For ku = V = T c /T the expression for — E j%'fP is, 
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the Poisson summation formula reads (/i = 2irm), 
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Noting that 



x 2 ln(l - e~ vx )dx = -C(4), 



integration by parts yields (||16|] p. 584), 
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Substracting F ^ T ^ gives Equation 20. 
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